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We apply the DMRG method to the 2 dimensional delta function potential which is a simple 
quantum mechanical model with asymptotic freedom and formation of bound states. The system 
block and the environment block of the DMRG contain the low energy and high energy degrees of 
freedom, respectively. The ground state energy and the lowest excited states are obtained with very 
high accuracy. We compare the DMRG method with the Similarity RG method and propose its 
generalization to field theoretical models in high energy physics. 
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A hallmark of an asymptotically free theory such as 
QCD is that it contains many degrees of freedom, with 
very different energy scales, which are coupled by the in- 
teraction Hamiltonian. Perturbative methods are valid 
for short distance physics but they fail for small momen- 
tum transfers or for energy scales where the bound states 
are formed. The existence of multiple energy scales sug- 
gests that the Renormalization Group approach is the 
correct strategy to attack these non perturbative prob- 
lems. In recent years there has been several proposals 
to extract effective low energy Hamiltonians using RG 
methods. Of particular interest is the light-front Hamil- 
tonian approach advocated in references which uses 
a similarity RG method (SRG) HQ. In this method the 
RG flow is given by an unitary transformation which di- 
agonalizes the Hamiltonian by succesive elimination of 
the off-diagonal matrix elements. The SRG-cutoff can 
be seen as the width of the band which contains the non 
vanishing off-diagonal matrix elements of the Hamilto- 
nian. At the end of the SRG flow the width is zero and 
the corresponding Hamiltonian contains in its diagonal 
all the eigenvalues of the original one. 

In this letter we shall propose an alternative RG ap- 
proach to study asymptotically free models using the 
Density Matrix Renormalization Group (DMRG). We 
shall also show the relations and differences between the 
DMRG and the SRG methods. The DMRG was proposed 
by White in 1992 to solve the problems of the old real 
space RG methods encountered in the 70 's, which led in 
those days to their abandon in favour of Montecarlo tech- 
niques H . The DMRG has by now become a standard 
numerical RG method applied to many body problems 
in Condensed Matter and other branches of Physics ( see 
references HQ for reviews). It is thus challenging to test 
how the DMRG handles the subtle dynamics of asymp- 
totically free theories. To our knowledge this is the first 
paper devoted to the subject. For this reason we have 
choosen as a theoretical lab a simple model possesing the 
essential properties of asymptotic freedom and formation 



of bound states, which are shared by realistic theories like 
QCD. 

The natural candidate for such a simple model is pro- 
vided by a 2d quantum mechanical particle subject to 
a delta function potential ||. The solution of the 2d 
delta function Schrodinger equation requires a regular- 
ization and renormalization schemes as in an ordinary 
quantum field theory. We shall use for our purposes the 
lattice regularization introduced by Glazek and Wilson 
in their SRG study of the problem ||,|l0| . These authors 
formulated the problem in momentum space where the 
states are labelled by an integer n that ranges between 
an infrared cutoff M and an ultraviolet cutoff N (i.e. 
M < n < N) . The kinetic energy E n of the state n in- 
creases exponentially as E n = b 2n , where b is an arbitrary 
constant greater than one. For numerical computations 
we shall take the value b = \/2 as in references [p|,^0[. 
The interaction Hamiltonian between the states n and m 
is given by —gy/E n E m , where g is the coupling constant 
of the problem. The discrete lattice Hamiltonian H is 
defined by the matrix elements 



L2n 



gb n+m , M<n,m<N 



(1) 



An overall shift of the levels by a constant term, i.e. n — > 
n + uq, implies that H nm scales with the factor b 2n °. 
This is a discrete version of scale invariance, which is 
broken by the infrared and ultraviolets cutoffs M and N. 
The latter symmetry implies that all the scales contribute 
to the observables, which makes very hard an accurate 
determination of their value by methods other than the 
exact one. 

The first step in the DMRG method is the partition of 
the system in two pieces called the system block and the 
environment block H. In our case we shall choose the 
system block Bf to be given by the low energy levels n 
which lie between the infrared cutoff M and the scale I 
(i.e. M < n < £), while the environment block Bf will 
contain the high energy levels n between the ultraviolet 
cutoff N and the scale t (i.e. I < n < N). The whole 
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system, with energy levels ranging from M to N, is ob- 
tained as the "superblock" Bf • oBf +3 , where • and o 
are the n = £ + 1 and n — 1 + 2 energy levels respectively 
(see fig.l). 
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FIG. 1. Superblock decomposition of the energy scales. 

The parameter £ varies from M to N — 3 and it labels 
the DMRG flow. Let us supose we want to find the GS 
of the whole system. We shall choose a trial GS wave 
function iptin) as follows, 




M <n<£ 
n = £ + l 
n = £ + 2 
£ + 3<n<N 



(2) 



where Li ( resp. Ri+z) is a normalized vector which 
describes the contribution of the low ( resp. high) energy 
block Bf ( resp. BEo) to the GS of the superblock B\ • 



oBf? +3 . The ansatz (|2j) is the momentum space version of 
the real space DMRG applied by White to study a free 
particle in a box . Our approach is close in spirit to 
the momentum space DMRG method proposed by Xiang 
p| . The energy of the state (|J) can be conveniently 
written as 



(1>t\H\tl> t ) = (a\H SB (t)\a) 



(3) 



where |a) is the vector (ai, <Z2, 03, 04) and the superblock 
Hamiltonian Hsb(£) is the 4x4 matrix given by 



( h L h Lm h Lo h LH \ 
hi, h, h, hn» 

\ hiH fin, hff hn J 



(4) 



The variational nature of the construction gives an upper 
bound of the exact GS energy. If the vectors Lg and Re+3 
coincide with the low energy and high energy pieces of 
the exact GS wave function then the DMRG algorithm 
presented so far would reproduce the exact result. Of 
course this is not in general the case but nevertheless, 
one can actually use the DMRG algorithm to improve 
in succesive steps the GS energy. The idea is to apply 
a continuity argument. Suppose we shift the scale £ to 
the next high energy level, say £ + 1. Then the new low 
energy vector L' e+1 will be related to the previous one Li 
by the equation 



whose entries read 

h L = {L t \H\L t ), h„ = (R e+3 \H\R e+3 ) 

h, = Ht + i,t + i, h = h^ + 2i+2 

hh» = J2 n =M Hni+\LfX n )i h Lo = J2n=M Hn,i+2Lt(n) 
h-Hm = En^+s^^+^wWi ^lh = (Le\H\Ri +3 ) 
hfio = Yl n =e+3 HnA+2Ri+z{n) h, = He +2 J+3 

(5) 

where H niTn are the matrix elements given in eq. No- 
tice that eq.(||) takes the form of an eigenvalue problem 
in a reduced vector space with only 4 degrees of free- 
dom. The GS of the superblock can be found by looking 
for the lowest eigenvalue E\{£) of the 4x4 matrix Hsb- 



M <n< 
n = £+l 



(G) 



where (a[, a' 2 ) — (oj., 0,2)/ \/a\ + a\ is the normalized two 
component vector obtained by the projection of the low- 
est eigenvalue of Hsb(^) into the block Bf». Similarly 
the energy h' L (£+l) associated to the latter block is given 
by 



ti L (£ + l) = (a' 1 ,a' 2 ) 



h L (£) h L .(£) 
h L .{£) h.{£) 



(7) 



The data L' e+1 and h' L (£+ 1) fully characterize the new 
block B'^ rl which can be regarded as the renormalization 
of the block B\». The next step is to construct the su- 
perblock Bf +1 »oBf +i which by the same techniques leads 
to the construction of a new block B' e ^_ 2 and so on. This 
procedure is iterated until the scale £ = N — 3, where one 
reverses the DMRG steps in order to update the high en- 
ergy blocks Bf using the low energy blocks built in the 
previous steps. After a few sweeps from low to high en- 
ergy and viceversa the lowest eigenvalue of the superblock 
Hamiltonian (Q) converges to a fix value which gives the 
DMRG estimation of the GS energy. To start out the 
process one has to grow up the system to its actual size. 
This can be done by considering superblocks of the form 
■Bm+p • °B§_ p where p = 0, . . . , (N - M- 3)/2. The last 
value of p yields a system containing all the scales from 
M to N. The low and high energy blocks constructed 
in the warm up are the starting point for the sweeping 
procedure explained above (see for details). The 
previous algorithm has been generalized in reference |l2| | 
to find out not only the GS but the low lying excitations 
as well. 

Let us now present our DMRG results for the case 
considered in references |9|JlO|], where M = —2l,N = 16 
and g = 0.06060600032108866. The latter value of g is 
choosen in such a way that the exact ground state energy 
of (Q) is given exactly by —1. The DMRG algorithm 
presented above gives the exact ground state energy with 
an error of 10 -14 (see table 1). Using the extension of 
the DMRG proposed in JTJ we have also computed the 
GS and the lowest 3 excited states of the Hamiltonian 
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(|). In table 1 we compare our DMRG results with the 
exact ones in terms of the relative deviation 



6E n = 



E n (DMRG) - E n (exact) 
E n (exact) 



(8) 



n 


1 


2 


3 


4 


5E n 


7 x lCT lb 


1.04 x 1CT Y 


3.36 x lCT b 


1.41 x lCT b 



Table 1. Relative error 8E n of the four lowest 
eigenstates of the Hamiltonian ([j]) . 

As shown in table 1 the accuracy of the excited states 
energies is lower than that of the GS. This feature is 
peculiar to the delta function Hamiltonian and it does 
not arise for the quantum mechanical models studied in 
[Il2j . There are several reasons for the very high accuracy 
of the DMRG applied to the Hamiltonian (jl]): i) the 
DMRG gives a variational upper bound to the exact GS 
energy which is usually improved in every DMRG step; 
ii) all the matrix elements of the whole Hamiltonian are 
used many times to feedback the superblock so that no 
information is lost; hi) the DMRG method focus on the 
determination of the GS and the low lying states. 

In fig. 2 we plot the DMRG wave function which after 
the third sweep is indistinguishable from the exact one. 
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FIG. 2. DMRG v.s. exact wave function. 

It is interesting to investigate the nature of the DMRG 
flow as compared with the one of the similarity RG 
method. In the SRG the effective Hamiltonian H(s) 
evolves as a function of s according to the Wegner equa- 
tion i], 



ds 



= [[H d ( s ),H(s)},H(s)} 



(9) 



where Hd(s) is the diagonal part of H(s). The initial 
condition of eq.@ is H(0) — H, where H is the origi- 
nal Hamiltonian of the problem. The parameter s ranges 
from to oo and it can be identified with the inverse 
square of the energy width A, i.e. s = 1/A 2 . Eq. (|) 



implies that H(s) is related to H by an unitary transfor- 
mation and therefore they share their eigenvalues. When 
s increases, the off diagonal matrix elements of H(s), lo- 
cated at distances greater than the width A = 1/ -^i, be- 
come very small. When s — oo the effective Hamiltonian 
-ff(oo) is diagonal and all its entries coincide with the 
eigenvalues of H . The numerical integration of eq.(^|) re- 
quires of course to follow the evolution of all the entries 
of H. One would like instead to project the effective 
Hamiltonians to smaller ( "window" ) matrices in order to 
reproduce the bound state eigenvalue jl0| . In a sense, the 
superblock 4x4 matrices (^) resemble the window matri- 
ces of ref . |l(| . Motivated by the SRG ideas |l(| , we have 
studied the RG flow of the eigenvalues Ei(£)(i = 1, . . . ,4) 
of the superblock Hamiltonian ([|). In fig. 3 we plot the 
lowest eigenvalue E\, together with the remaining ones 
scaled down by a factor b 21 . We can clearly see from fig. 3 
that Ei stays constant through all the DMRG steps while 
Ei(i = 2,3,4) vary with the energy scale b 21 with some 
deviations depending on the energy region. The plateaus 
correspond to low energy regions while the oscillations 
and bumps occur for intermediate and high energies. To 
a first order approximation, which is almost exact for the 
plateaus, the superblock Hamiltonian (|]) can be written 
as 



H SB {t) = O i 



( E 1 



V o 





E' 2 b 





2' 






E' 3 b 21 








E' A b 2e 



01 



(10) 



where Oi is a unitary matrix. Using eq.( 10) one can show 
that the superblock Hamiltonians satisfy the following 
second order recursion relation, 



b+b- 



H SB (i) = 

b- x U t H SB {l + l)U\ + bUj^Hseii - l)U t . 



(11) 



where Ui = Of O\ +1 . The continuum limit of cq.(|Tl]) 
gives the flow equations 



Hi 



dH SB 
dl 



[V,H. 



SB\ 



(12) 



where r\ = ^jfO\. Eq.([T^) is a second order diferential 
equation which is to be compared with the first order 
equation eq. ([)]) . The DMRG flow is a sort of similarity 
transformation with some eigenvalues running with the 
scale. Using the standard RG terminology the lowest 
eigenvalue E\ can be associated with a marginal operator 
while the eigenvalues Ei for i = 2,3,4 are associated 
with infrared irrelevant operators which vanish at the 
fixed point Hamiltonian Hgg(£ = M). Indeed all the 
entries of Hsb{@ — M) are very small except for the 
entry h B — —0.999 whose value is close to the bound 
state energy. These results suggest that the exactness 
of the DMRG method is due to a careful treatment of 
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the irrelevants operators, which in other RG methods 
are difficult to control in general. 

From a conceptual point of view the DMRG offers a 
new way of thinking about cutoffs and RG flows in high 
energy physics. Traditional cutoffs remove high energy 
states while the lowering of the cutoff produces effec- 
tive operators for lower energies [i3fl . In the Lagrangian 
formulation this strategy can be implemented perturba- 
tively without much difficulty However in the Hamilto- 
nian formulation it gives rise to small denominators prob- 
lems involving energy differences between the states kept 
and the states truncated in the RG process p^,p[. This 
latter problems do not arise in the DMRG truncation 
for it uses a non perturbative self-consistent algorithm to 
find the best choice of the effective Hilbert spaces and 
Hamiltonians. 

The next step in the application of the DMRG to high 
energy physics is of course to consider field theoretical 
models with asymptotic freedom and bound states. 
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3. Plot of the rescaled superblock eigenvalues Ei/b 21 
3, 4) as a function of the DMRG step L from sweep 2 



The most appropriate formalism for this application of 
the DMRG is known as Discrete Light-Cone Quantiza- 
tion in momentum space (DLCQ) [ fl5| |, ]l6[ |. In the DLCQ 
approach the Hilbert space is finite dimensional and the 
light-front Hamiltonian Hlc acting on it is similar to 
that of a many-body Hamiltonian in condensed matter 
p6[ . The search for bound states amounts to solving the 
Schroedinger equation 



H 



LC\ 



= M 2 



(13) 



where M 2 is the mass of the bound state. Thus, one 
can apply to ( |l3|) standard diagonalization techniques 
such as Lanczos. The DMRG method allows us to study 
larger Hilbert spaces than those achieved with the Lanc- 
zos method. This is needed in order to recover the con- 
tinuum limit of Hlc- 

The key to make the clever truncation of states in the 
RG process is given by the density matrix of the blocks. 



This is a systematic procedure with no wild guessing 
about the wave function. A given block will contain the 
most representative multiparticle and momentum states 
in order to reconstruct the lowest lying bound states. On 
the other hand, the DMRG is a numerical method, un- 
like the more analytical SRG method |l7|, and does not 
need perturbative inputs. The basic requirement for the 
DMRG method to work is a discretized Hamiltonian act- 
ing on finite dimensional Hilbert spaces. As a first test of 
this program we have solved eq. ( [j~3| ) for the Positronium 
state of the massive and massless Schwinger model in the 
one-fermion sector with the DMRG. These results will be 
reported elsewhere. 

Thefore, we expect that the main ideas presented in 
this letter can be generalized to the DLCQ Hamiltoni- 
ans. Specifically, the breaking of the system into low 
energy and high energy blocks which are constantly up- 
dated through the DMRG process. On the other hand 
the DMRG method combined with the DLCQ approach 
does not have the sign problems that emerge in the Mon- 
tecarlo methods used in Lattice Gauge Theories [jl(|. In 
summary, we believe that the power shown in condensed 
matter systems by the DMRG method is worthwhile to 
be translated into Particle Physics. 
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